Response of Unruh-DeWitt detector with time-dependent acceleration 
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It is well known that a detector, coupled linearly to a quantum field and accelerating through 
the inertial vacuum with a constant acceleration g, will behave as though it is immersed in a 
radiation field with temperature T = (g/2n). We study a generalization of this result for detectors 
moving with a time-dependent acceleration g(r) along a given direction. After defining the rate 
of excitation of the detector appropriately, we evaluate this rate for time-dependent acceleration, 
g(r), to linear order in the parameter n = gig 2 ■ In this case, we have three length scales in the 
problem: g —1 , (g/g)^ 1 and where uj is the energy difference between the two levels of the 
detector at which the spectrum is probed. We show that: (a) When uj^ 1 <C g 1 <C (s/sO -1 ! the 
rate of transition of the detector corresponds to a slowly varying temperature T(r) = g(T)/2n, as 
one would have expected, (b) However, when g' 1 <C uj^ 1 <C (<?/(?) _1 , we find that the spectrum 
is modified even at the order 0(n). This is counter- intuitive because, in this case, the relevant 
frequency does not probe the rate of change of the acceleration since (g/g) <C u> and we certainly 
do not have deviation from the thermal spectrum when g = 0. This result shows that there is a 
subtle discontinuity in the behaviour of detectors with g = and g/g 2 being arbitrarily small. We 
corroborate this result by evaluating the detector response for a particular trajectory which admits 
an analytic expression for the poles of the Wightman function. 

PACS numbers: 04.62.+v,04.60.-m 
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I. INTRODUCTION 

One of the key results which emerge from the study 
of quantum field theory in non-inertial coordinate sys- 
tems (and curved spacetime though we will not consider 
it in this paper) is that both the particle content of the 
quantum states, as well as the pattern of vacuum fluctu- 
ations, are not generally covariant. This can be explicitly 
demonstrated by studying the response of detectors lin- 
early coupled to the quantum field (usually called Unruh- 
DeWitt detectors) in different states of motion ([1, 2]; see 
[3] for a review). The probability that a detector trav- 
eling along the trajectory XJ(r) will get excited can be 
expressed as an integral over the Wightman function of 
the field in the form 
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dr 2 / dn exp (-iurn) G + [Xt(T 2 ),Xi( Tl )] 
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dt / du exp (— iuju) t] 



(i) 



where uj is the energy difference between the two levels 
and the coordinates (u, t) are introduced through the def- 
initions: u = T2— T\, It = T2+T1, and we have absorbed a 
factor dependent upon the internal details of the detector 
in the definition of V . 

When the trajectory X\(t) is along the integral curve 
of a time-like Killing vector field in flat spacetime (we will 
call such trajectories 'stationary'), the Wightman func- 
tion G + [u, t] will only depend on the time difference u so 
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that G + [u,t] = G[u]. Then the second integral over t in 
Eq. (1) will lead to a divergent result. This is handled 
by the usual procedure of time dependent perturbation 
theory which involves ignoring the integral over t and 
interpreting the rest of the result as providing the rate 
of transitions between the two levels. For the stationary 
trajectories, this rate will be a constant. 

It is also worth mentioning at this point that the rate 
so defined is a real number. This is easily seen by noting 
that, the imaginary part of the rate will be given by: 



Im[£>] = 



du < exp (— iuiu) G + [u, t] (2) 



exp(iuju) (G + )*[u 7 t] 



We now note that the Hadamard function satisfies: 
(G+)*[P,Q] = G+[Q,P] for any two points (P,Q), here 
characterised by (n, T2). Further, it follows from the def- 
initions of u and t that interchanging the points amounts 
to u — } — u and t ->• t, so that (G+)*[u,i] = G + {-u,t]. 
It is then easy to see that lm[P] — 0. 

To avoid possible confusion, we must also point out 
that our definition for the response function differs from 
certain other definitions found in literature (often mo- 
tivated by arguments of causality etc.). Essentially, the 
difference lies in the choice of the "time" variable with re- 
spect to which the rate is defined (which, in our case, is t). 
Similarly, our choice of regularization scheme is also dif- 
ferent from some other choices found in literature. Hence, 
the result we shall obtain can not be directly compared 
with other results based on a different choice of definition 
for the rate or regularization scheme. We shall have more 
to say on this in the last section. 
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Since flat spacetime admits ten independent Killing 
vector fields, one can construct several linear combina- 
tions of these Killing vectors which will be time-like in 
parts of the spacetime. The response of detectors on 
these trajectories have been studied extensively in the 
literature (see, e.g., Ref. [4-6]). In a generic situation, 
the detector will respond to the pattern of vacuum fluc- 
tuations which can coincide with the particle content of 
the quantum field determined by Bogoluibov coefficients 
in specific cases but not always. ( e.g., for a circular 
trajectory, the particle detector "clicks" but the number 
of particles calculated using the Bogoluibov coefficients 
turns out to be zero. Many of these conceptual issues 
have been discussed and clarified in the literature [6-8].) 

The uniformly accelerated trajectory corresponds to 
the integral curve of the Killing vector field corresponding 
to the Lorentz boost along the direction of the accelera- 
tion g. In this case [9] we obviously have G + [u,t] = G[u}. 
In this particular case, the pattern of vacuum fluctuations 
match with the particle content of the quantum state 
and the rate of excitation of the detector will correspond 
to a thermal spectrum of particles with a temperature 
T = g/2n. This is of particular importance because it 
allows us to associate a temperature with the Rindlcr 
horizon with obvious implications for black hole physics. 

Unfortunately, a detector which is uniformly acceler- 
ated from r = — oo to r = +oo is not physically real- 
izable. The question arises as to what happens in the 
case of more realistic detectors. One possible way of ad- 
dressing this question is to keep the coupling to the field 
switched on only for a finite interval of time (see, e.g, 
Ref. [10]). But this introduces transients and one needs 
to handle them with care. It also does not seem very nat- 
ural to switch off the coupling in this manner. A more 
obvious and physically interesting way of attacking the 
problem would be to study the response of a detector 
moving along a given direction with a time-dependent 
acceleration g(r) which is what we will do in this paper. 

There are three further motivations for taking up this 
study which are somewhat indirect. 

First, we know that there is a direct correspondence 
between the detector response in a uniformly acceler- 
ated trajectory and the phenomena which takes place in 
the Hartle-Hawking vacuum state around a black hole. 
By extending this analogy, we would expect a sub-class 
of time-dependent accelerations — especially those g(r) 
which vanish at early times and become constant at late 
times — to correspond to the phenomena which takes 
place in a collapsing black hole scenario in the Unruh 
vacuum state. (For preliminary discussions along these 
lines, see section 5.1 of Ref. [11].) This would be inter- 
esting to study. 

Second, there has been considerable amount of work in 
recent years which attempts to interpret the field equa- 
tions of gravity as a thermodynamic identity. This body 
of work [12] uses the concept of local Rindler observers 
that corresponds to trajectories which, in the local iner- 
tial frames around any given event, will be a hyperbola. 



While one expects such a local concept to be valid as a 
first approximation, it is important to verify it explic- 
itly (and indeed our results in this paper will justify this 
notion and make it sharper). 

Finally, this subject has thrown up fair number of sur- 
prises and subtleties in the past and one cannot take it 
for granted that intuitively obvious results will arise when 
we rigorously analyse the case of, say, a slowly varying 
acceleration! It requires explicit verification. Our naive 
expectation will be that, for sufficiently slowly varying 
acceleration (with {g/g 2 ) -C 1) one would expect the de- 
tection rate to correspond to a time dependent tempera- 
ture T(t) cx g(r). At the same time, one will not expect 
such a result to hold for all frequencies of the thermal 
spectrum. There is, in fact, a good reason to expect some 
modification due the presence of (local acceleration) hori- 
zon. This sets a length scale g -1 in the problem, which 
can be compared with the length scale probed by a par- 
ticular mode, tu^ 1 . Of course, we know that the spectrum 
is Planckian for all values of when g is constant; it is 
therefore interesting to see whether a varying g makes any 
difference. As we shall show, one does get low frequency 
(g ui <c 1) modifications when g is non-vanishing even 
when g/g <C lo, which is a surprising result. 

In Sec. II, we describe the setup appropriate for calcu- 
lating the response function. In Sec. Ill, we evaluate the 
Unruh-DcWitt detector response for time-dependent ac- 
celeration, <?(t), to linear order in the parameter rj = 
g/g 2 . We find that, to this order, the spectrum can 
indeed be approximated in the UV region (cj 3> g) 
by Planck spectrum with time-dependent temperature, 
T = g(r) /2tt. However, the spectrum is modified even at 
0(rj) for to <c g. In Sec. IV, we corroborate this result 
by evaluating the detector response for a particular tra- 
jectory which admits an analytic expression for the poles 
(under a particular approximation). Finally, we conclude 
with few relevant comments. We use the metric signature 
(-,+,+,+). 



II. DETECTOR RESPONSE: BACKGROUND 

The trajectory of an observer moving with a time- 
dependent acceleration, <?(t), with r being the proper 
time, is given by 

T*(t) = J da coshx(a) 

X*(t) = J da sinhx(a) (3) 

where 

^~3(t); X (r) = £ da g(a) (4) 

and (A ,^ 1 ) = (T, X) arc inertial coordinates. (We 
have chosen x(0) = to obtain the integral form.) The 
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local coordinates of the observer, (r, x), can be con- 
structed easily; these are given by (see, for e.g, Equation 
(73) of Ref. [11]) 



T(t) = / da [1 + g(a)x] coshx(a) 

X(t) = [ da [l + g(a)x] sinhx(a) (5) 



In the local coordinates, the observer is always located at 
x = 0. We shall assume T2 > T\ without loss of generality 
The probability of transition for the detector is given 
by Eq. (1) which is valid for any trajectory. In general, for 
an arbitrary g{r), there is no time translational symme- 
try and G + will depend on both u and t. Following the 
procedure adapted for stationary trajectories, we shall 
again define the transition rate by ignoring the integral 
over t. But now this rate will be time-dependent, due 
to the t— dependence of G + , which, of course, is to be 
expected. So, we shall define the transition rate to be 



V 



duexp (— iuiu) G + (u,t) 



(6) 



The Wightman function is given by 

« + (M) = JLi 

£ 2 (1,2) = -[T.fo)-T*(Ti)] 2 + [X.to)-X*(Ti)] 2 

(7) 

with an ie prescription which is implicit in the difference 
of the time coordinates. Substituting Eqs. (3), the ex- 
pression for £ 2 can be written in the following convenient 
form; 



-U I- 



where 



pT2 

I±= da exp±x(a) 
J n 



(8) 



(9) 



We see that the detector response is essentially deter- 
mined by the poles of £ 2 . For the constant acceleration 
case, the poles of I + and I_ coincide, so that we have an 
infinity of second order poles, the residues at which gives 
the well known thermal response function. Our task, 
therefore, reduces to identifying the poles of I± and eval- 
uating the integral in Eq. (9). We shall now turn to this 
task. 



III. DETECTOR RESPONSE FOR g/g 2 < 1 

As one can easily see, it is impossible to determine the 
structure of the poles for a general g(a). Hence, we shall 
attack this problem in two steps. First, in this section, 
we will consider a slowly varying acceleration and obtain 



the detector response. In the next section, we shall work 
out the response for a specific form of g{r). 

Consider a general g(r), which varies slowly compared 
to its value go at r = which can be chosen to be an 
arbitrary instant of proper time. We shall now expand 
g(r) in a Taylor series retaining the lowest order terms: 



9(r) = 
X(r) = 

Therefore, we have 



g r + O(g T ) 
1 



So - 

9or + -g T z + O{g T A ) 



9ot 



1 + r - 



exp±x(r) = exp (±g f) 



1±-7 ? o(5ot) 2 +0(t? 2 ) 



(10) 



(11) 



where we have defined tjq = go /g^ , and we shall do subse- 
quent calculations keeping terms up to O(?7o). A trajec- 
tory is, of course, not completely specified by r)Q. For our 
result to remain valid, the contribution of higher deriva- 
tives of acceleration must be ignorable compared to g. 
Although a restriction, this condition will almost always 
be fulfilled in physically relevant cases, when there is only 
one small parameter in the problem. If not we will get the 
same result when all the corresponding higher derivatives 
of the acceleration are small. 

We now proceed to analyze the pole structure of £ 2 
to determine the detector response. Evaluation of I± 
involves trivial integrations; we obtain, 



I± 



1 

go 



l d 2 



where we have defined 



Q± (£1.62; ") 



Q± (6, 6; ") 



d£ exp±a£ 



(12) 



(13) 



with £ 1(2 ) 
1 



: 50^1(2) 

gl 



Q+Q- 



Therefore, we obtain 

, m 1 



2 Q+Q- 



(Q+Ql-Q-Q'i) 



sign implies that we have ignored 0(r/g 



(14) 
terms, 



where 

as we should for consistency. (The prime stands for d/da, 
with a being set to unity in the end.) 

The zeroth order term is just the constant accelera- 
tion Rindlcr contribution. We shall now analyze the pole 
structure of the second term. This term can be further 
simplified using expressions for Q±. Specifically, the term 
in the round brackets in Eq. (14) can be written, as 



Q+Q- - Q-Q+ =2AQ-- 2B 



(15) 



where 



A 
B 



Kcxp+^--[£ 2 cxp+£]^ 

[£exp-e]! + i[£ 2 cxp-^; 



(16) 
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with the obvious notation: 



■K60 



•]&)■ 



(17) 



It can be shown that A and B are both finite at the 
poles. The second term in Eq. (14) therefore has cubic 
order poles determined by zeros of Q±. Substituting the 
above expressions into Eq. (14), we obtain 



1 



0+0- 



0+0- 



A_ 

Q+ 



B 



(18) 



The first term is the standard Rindler contribution, and 
it is well known that this term gives a second order pole 
at Uk = —iftok with k > 0, and /3q = 2tt / g$. From 
here onwards, we shall denote derivatives with respect to 
u by an overdot. We note that Q+Q-jg^ = 1 at the 
poles, which is the standard result for Rindler and can 
be easily verified by explicit computation (usually, one 
uses the well known infinite image sum representation of 
(sinha;) -2 to obtain the same result). The pole structure 
is now determined by 



Q+Q-Q± = Q+Q-Q± (u - u k f + 0((u- u k f) 



(19) 



where the u derivatives are to be evaluated at u k . To 
compute the residues, we need to evaluate second deriva- 
tives with respect to u of the functions Aexp —iuju and 
i3exp— iuu, at u = Uk- This is straightforward and we 
relegate the details to Appendix A. The calculations 
are enormously simplified by noting that, as t — > —t, 
B — >• A, so that we need to consider only terms odd in t 
in Eq. (18); the remaining terms (which would otherwise 
be tedious to evaluate), cancel. 

The transition rate of the detector is given by Eq. (6) 



duexp(— iuju) G (u,t) 
exp (— iuu) 



1 

'4^2 



dv 



(20) 



Substituting the residues at the poles, calculated in Ap- 
pendix A, we obtain 



1 °° 

V = — uj exp —ftoLuk 



k=l 



TT 

(3 uj 



k exp —ftpuik (21) 



fe=i 



which is correct to 0(770). Hence, we see that the resul- 
tant spectrum will not be thermal at all frequencies even 
to order O(r?o), due to the second term in the square 
brackets, which becomes significant at low frequencies. 
[A similar result was arrived at recently in [14] in a dif- 
ferent physical context.] 



In the UV region (i.e., Pquj 3> 1), we get 
1 ^ 



V 



2tt ^ 

k=l 



u> exp —(3 uik 



Noting that 5/3 
V 



1 

2^ 



+ (i]ot)uj 2 k exp — /3owfc 

k=l 

= —2irr) Q t, this can be written as 
d 



(22) 



1 + 5(3 



d(3 



uj exp —j3ujk 



1 

27 



1 



p=p \exp/3w 



(23) 



Therefore, to 0(770), we have 

1 UJ 



V = 



2tt exp [f3(t)uj] 



1 



— uj exp 
2tt 1 



■WW (24) 



where /3(t) = 2w/g(t). This result is intuitively under- 
standable; at sufficiently high frequencies, we just recover 
the usual result with g replaced by g(r) when the accel- 
eration varies with time. 

However, note that it is valid only in the UV region; 
our result also shows that the spectrum will be modified 
for (3 uj <C 1. In fact, the second sum in Eq. (21) is easily 
evaluated, and we obtain 



P = I P [g ] + vtco 2 



i-C- 

s 



exp[s] 



[cxp[s] - 1]' 



(25) 



where, for convenience, we have defined s = /3qui, 
and Ip[go] represents Planck spectrum at temperature 
go/(2n). As stated above, for s ^> 1, the second term 
in square brackets can be neglected and the remaining 
terms combine to give V ~ Ip[g(t)}. 

We want to analyze the s < 1 case a bit further, to 
highlight a counter-intuitive fact. In this limit, we obtain, 



V 



ip\go}-v tuj (- 



ip[go] - vtu 2 (-0 



In the same limit, ip[go] ~ u/(2tts) = 
we can rewrite the above expression as 



(26) 

C (say), so that 



C 



1 



(2. 3 )4- 



= C [1 - (2ir 3 )ba 2 ojt\ (27) 

where b and a are the dimensionless quantities, 

a = 1/s = g a /27ruj; b = g /(g uj) (28) 

Evidently, there are two possibilities, depending on 
whether b is greater than or less than one. When 



s <C r/ <C 1, or 
1 < b < a 



(29) 
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we see that the frequency is probing the change in accel- 
eration because cu^ 1 3> (g/g) ■ So we certainly expect 
the spectrum to be distorted and this is what happens 
and the result is understandable. However, when 



?7 < s « 1, or 
b < 1 < a 



(30) 



we see that g/g <C w and so these frequencies are not 
probing the change in acceleration at all. Therfore, one 
would have expected to recover the results of constant 
acceleration in which case there are no distortions from 
the thermal spectrum at any frequency. But we see that, 
in this case, we can still have ba 2 <C 1 and produce a 
distortion of thermal spectrum at low frequencies. 

Before proceeding further, we must highlight an impor- 
tant assumption that has gone into the derivation. For 
calculating the residues, we have first expanded the inte- 
grand and then evaluated the residues. The true expres- 
sion is, of course, to be obtained by first doing the contour 
integral and then expanding in r\. We have assumed that 
the two steps, Taylor expansion and integration, com- 
mute. Our result will be invalidated for functions g{f) 
which fail to satisfy this criterion. Moreover, we have also 
assumed that the u integration goes all the way from — 00 
to +00, while doing a Taylor series in t. It is important 
to understand better whether such an approximation is 
valid, and, if not, what difference will it make to the re- 
sult. In particular, the low frequency modification we 
obtain may be an artifact of such a truncation, and this 
caveat must always be kept in mind. 



IV. DETECTOR RESPONSE FOR A SPECIFIC 
TRAJECTORY 

We shall now study the response function correspond- 
ing to a particular detector trajectory determined by 



g{r) 



go 



l + eg \T\ 



(31) 



where e is a small, dimcnsionlcss parameter. The re- 
sponse can now be evaluated in a straightforward man- 
ner. We expect a splitting of the quadratic poles at O(e) 
from their constant acceleration values, so that we es- 
sentially have a couple of first order poles separated in- 
finitesimally. For the above g(r), we have 



X(r) = (l/e)ln(l + e. 9o |T|)sgn(r). 



(32) 



Because of the dependence on |r|, we need to consider 
the cases (i) < t\ < t 2 , (ii) t\ < r 2 < and (iii) 
ti < 0, r 2 > separately while evaluating I± (we take 
T2 > T\ without loss of generality). While the first two 
cases admit analytic expressions for the poles, the same is 
not true of (iii). We shall evaluate the response function 
for the first two cases, and comment on the possible effect 
of (iii) later on. 



The integrals involved in I± are trivial; for clarity, 
we refer to values of I± for case (ii) as I±, and those 
for case (i) simply as I±. Then, it is easy to see that 
I±(t\,T2) = T2,— Ti). So, we can obtain I± from 
I± simply by changing t to —t (see definitions of u and t 
above), or, what is the same thing, by replacing t with \t\ 
in the expression for case (i). With this understanding, 
we simply write t rather than \t\ in the expressions below. 
We also define £1(2) = 1 + e<?oTi(2) and 77 = £2/^1 • Then, 



: ±(l/e) c ±(l/e) 



(33) 



The poles (i.e., the zeros of I±) are determined by: (1 ± 
e _1 )ln?/ = 2nik. Now rewrite £1 and £2 in terms of u 
and t, to obtain 



2i 



-(e-i+got) tanf^ 

go V 2 



(34) 



where xt = 2nk/(l ± e^ 1 )- As a check, note that 
this reduces to the standard constant acceleration val- 
ues, ±2iirk/go for e = 0. 

Rest of the calculation involves standard residue cal- 
culus, and is quite lengthy Since the poles are now 
split at O(e), we need to evaluate the quantities R+ = 
I' + (u~^)I_(u~^) and i?_ = I'_(uu)I+(uu) for calculating 
the residues. The quantity R± x go/0- + £<7oO is given by 



1 + cxp 



±i- 



exp [ixt (1 T \)] 



ex P[«Xfc] 



(-e±l)(l + exp[i X ±)] 



(35) 



It is now a straightforward exercise to use Eqs. (34) and 
(35) and evaluate the response function. This turns out 
to be 



1 00 
2tt ^ v 



k=i 



No further calculations are required, since it easy to see 
that, with Sg = —eg 2 t + 0(e 2 ), the two terms above com- 
bine to give a Planck spectrum with temperature g(t) /2ir. 

Let us now turn to the contribution of poles which we 
have not accounted for. In the above calculation, we left 
out the contribution to the integral of the it-range where 
Ti,T2 have opposite signs. Unfortunately, this case does 
not admit analytic expressions for the poles. But, from 
the result in Sec. Ill, we expect that this contribution 
will be irrelevant at high frequencies. Apart from this, it 
is not possible to make any comments about this contri- 
bution. As already mentioned in the Introduction, this is 
typical of most of the calculations that attempt a rigorous 
evaluation [in particular, [13] discusses characteristics of 
detector response in curved spacetime] although the ex- 
plicit result we have obtained is very close to what one 
would have expected for a slowly changing acceleration. 
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V. CONCLUSIONS 

The result brings out another interesting fact associ- 
ated with the combined effect of presence of the horizon 
and varying acceleration and we shall discuss this briefly. 

We see that, for modes with co^ 1 <C the result es- 
sentially involves replacing the acceleration by its instan- 
taneous value so that T(r) = 5(r)/27r but the thermal 
spectrum gets distorted for uj^ 1 3> g -1 . At first sight, 
one would have thought that this is to be expected. We 
know that for an accelerated trajectory, g~ l gives the 
approximate location of the local horizon. (The exact lo- 
cation of the horizon will change with r, see Appendix B). 
On the other hand, a mode with frequency lo will probe 
a length scale ~ oj" 1 . Such a mode will be within the 
region 'outside the horizon' if u^ 1 <C g^ 1 , or f3uj 3> 1 but 
will probe the horizon scale and beyond if 3> g^ 1 ■ 
So one may think that it is natural for the spectrum to 
be distorted in the latter case. 

There is, however, a subtlety here. Our problem actu- 
ally has three length scales not just two: and 
quite crucially, (<?/<?) It would have been no surprise, 
if the spectral distortion arose for lu^ 1 ^> {g/g)^ 1 ; these 
are the frequencies which see the change in the accel- 
eration and there will be some distortion. But this is 
not what we found! Instead we find that in a typical 
situation with g~ x <C (g/g)~ l , with very slowly varying 
acceleration, the spectral distortions occur already when 
g~ l <C uj^ 1 <C (g/g)^ 1 . So whether g ^ or whether 
g = makes a difference to the spectrum even when the 
relevant frequency is not probing the time variation of 
the acceleration. Obviously, this effect does not exist in 
the case of g = 0; so we need to conclude that one cannot 
take the limit continuously for all frequencies. We believe 
this arises due to the changing distance to the horizon 
but only further investigations will nail down the precise 
reason. 

It is known that viewed from the inertial frame the 
final state of the field is an one-particle state so that 
the 'detection' is actually accompanied by an emission. 
It has been suggested [sec [7]] that it is better to think 
of the detector as radiating a Minkowski particle, rather 
than "detecting" anything. From this point of view, it 
would be interesting to see whether the response func- 
tion we have obtained has some simple interpretation, 
particularly the g term. 

Finally, as promised in the Introduction, we briefly dis- 
cuss the issue of comparing our results with other results 
in literature. In doing so, one must realize that our ex- 
pressions for the rate as well as the regularization are dif- 
ferent from the ones used in the literature. Our choices 
are based on the simple fact that it is closest to what one 
does in the standard, constant acceleration case. A dif- 
ferent choice of variables for defining the rate (and even 
a different regularization scheme) can alter the results 
since the pole structure will change. In this context, we 
would particularly like to mention the analysis presented 
in [13, 15] Our choice of regularization (with an it pre- 



scription on u) is actually similar to the one employed 
in [15] (see their Eq. (21)). However, the difference lies 
in the definition of the response rate itself. In particu- 
lar, the relevant function which gives transition rate at 
time t, for e.g., in Rcfs. [13] and [15], is G(t,t — u). It 
is not difficult to sec that this would have a completely 
different functional dependence on t and u as compared 
to our case, since the definition of t is manifestly differ- 
ent. So, effectively one is integrating completely different 
functions of u in the two cases; the results can not, there- 
fore, be directly compared as such. Additional physical 
criteria are needed to choose one definition over another 
(for e.g., in [13] and [15], the motivation is causality. In 
our case, the coordinates t and u corresponding to the 
two points on the trajectory actually correspond to the 
so called "radar coordinates" which are natural set of lo- 
cal coordinates assigned to nearby points connected to 
( T ii T 2) by light beams. By their very construction, these 
coordinates are non-local, and the transition rate must 
be interpreted keeping this in mind. However, further 
work is needed to make a precise connection.) 
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Appendix A: Evaluation of residues of Eq. (20) 

In this appendix, we outline the evaluation of residues 
at the cubic order poles of Eq. (20). We would only stress 
certain steps which are crucial to minimize the algebra 
in the otherwise elementary calculation. 

To begin with, define two new variables, a = gou/2 
and b = got, so that, £i = b — a and £2 = b + a. We then 
have 



A 



B 



2 e +h cosh< 



2 e cosh a 



a(l - b) + tanh a ( +b - -b 2 - -a 2 



a(l + b) + tanh a 



2 

-b--b 2 
2 



2 
1 , 



(Al) 



It is evident that B(a,b) = A(a, — b), as mentioned in 
the text. We essentially require [see Eq. (18)] A/Q+ and 
B/Q-. From their definition (13), we have 



Q± = .9o e ±b cosh a 



(A2) 



Note that Q+Q-fg 2 . — 1 at the poles, which, as was em- 
phasized in the text, is the standard result for Rindlcr. 

So, we have A/Q+ (a,b) = B/Q_ (a, -b). Therefore, 



7 



the contribution of the 0(770) term in Eq. (18) to V be- 
comes 

/ 1 \ f-2m\ 7? [ d 2 " 

(A3) 

plus a similar term for B, and a sum over all relevant fc's 
(such that Wfe's lie on the negative imaginary axis, with 
the contour closed in the lower- half complex it-plane). All 
that remains is to pick out the terms in A which are odd 
in b, calculate the second derivatives which are required, 
evaluate at u%, and then multiply by 2 for the contribu- 
tion of the B part. These are all straightforward, though 
lengthy, steps. The object of interest is [sec Eqs. Al, A2 
and A3], 

f = e~ b A = 2(cosha) [b (tanha — a)] + (terms even in&), 

(A4) 

and, at the poles, we obtain, 

/ = 0, / = — 2a/ c a 2 ,6 cosh cifc and / = —2akbcoshak- 

(A5) 

Putting everything together, we finally obtain Eq. (21). 
[Note that the second term in the square brackets in 
Eq. (21), which becomes significant in IR, arises from 
the / term above.] 

Appendix B: Local horizon for a trajectory with 
time dependent acceleration 

In the local coordinates based on a trajectory with 
time-dependent acceleration [see Eq. (5)], the metric be- 



comes (see for e.g, Rcf. [11]) 

ds 2 = - [1 + g{r)x} 2 dr 2 + dx 2 + dY 2 + dZ 2 (Bl) 

The equation for a null surface can be written as 
<&(t, x) = x — /(t) = 0. The function /(r) is deter- 
mined by the condition 9 a <I><9 a <I> = 0. Doing this leads 
to a differential equation for /(r); it's solution yields the 
following expression for the horizon location 

xu(t) = pexp-£,(r) J dy exp£(y) (B2) 

where 

£(r) = T / 9(x)dx (B3) 

For g{r) ps g + g r, we have 

Mr) = -g^ [1 - ?7(±1 + <?ot) + 0(r, 2 )] (B4) 

We can also invert this to write, to the same order of 
accuracy, 

g(r) ~ -Xh 1 [1 T v] (B5) 
The temperature T(t) — g(T)/2ir, associated with the 
detector response at 0(rj) and in the UV limit, can be 
cast in an interesting form by further noting that, v = 
±h = r/ + 0(?/ 2 ): 

T ( T ) - 4-^ - -7r~T~r~ m 

Zir Zirxii 1 ± v 
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